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Abstract. The effect of rotation of a star on the distri- 
bution of chemical species in radiative zones is discussed. 
Gravity darkening generates a large radiative force on 
heavy element ions which is directed toward the equato- 
rial plane. Taking iron as an example, it is shown that this 
force may produce drift velocities similar to, and larger 
than, the typical velocities of bulk motion due to merid- 
ional circulation. This potentially allows large chemical 
abundance inhomogeneities to build up across a meridian 
over the lifetime of the star - particularly near the equato- 
rial plane. This enhancement may be significantly reduced 
if the mass loss of the star is strongly metallicity depen- 
dent, in which case the mass-loss rate may be enhanced 
in the equatorial plane. 

Key words: stars: abundances, evolution, interiors, rota- 
tion - diffusion 



1. Introduction 

Abundance anomalies are observed in many different types 
of star - the photospheres of stars with effective tempera- 
tures hotter than ~ 6000K may have undergone chemical 
separation (Michaud 1970). It is generally accepted that 
this may be explained via the force on atoms by radiation: 
if this is larger than gravity, then the element will diffuse 
toward the surface, else it will settle. The force due to radi- 
ation pressure on elements in a star has been discussed at 
length by e.g. Michaud et al. (1976) Vauclair et al. (1978), 
Aleican & Artru (1990) and Gonzalez et al. (1995). Abun- 
dance anomalies have been discussed in a variety of chem- 
ically peculiar stars eg. Am and Fm stars (Charbonneau 
& Michaud 1991, Alecian 1996), A Bootis stars (Michaud 
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& Charland 1986) and magnetic Ap-Bp stars (Alecian & 
Vauclair, 1981, Michaud et al. 1981). When the star is also 
rotating, the surface abundance anomalies generated by 
radiation pressure may be destroyed by the action of tur- 
bulent motions within the star (Charbonneau 1992, Talon 
& Zahn 1997). These arise from shear motions created by 
the meridional motions which arise due to gravity dark- 
ening (see eg. von Zeipel 1924, Tassoul & Tassoul 1982, 
Zahn 1992) 

There is some observational evidence for latitudinal 
abundance variations of heavy elements in planetary nebu- 
lae (Balick et al. 1994). This is possibly explained by more 
highly processed material being ejected faster but at later 
times then the rest of the shell. However, the observations 
raise an interesting question: it is possible that there are 
stars which are chemically inhomogeneous in both radial 
and latitudinal directions? 

Here a non-magnetic mechanism which generates lat- 
itudinal abundance variations is suggested and investi- 
gated. It is found that a significant latitudinal radiative 
acceleration on heavy elements exists in rotating stars due 
to gravity darkening (von Zeipel 1924), which causes ionic 
species to diffuse toward the equatorial plane. Competing 
with this equatorward drift are the meridional circulation 
currents, and turbulence generated via shear motions. The 
aim of this paper is solely to assess whether any latitudinal 
drift of heavy elements may occur in rotating stars. 

The transport of chemicals is discussed in §2. In §3 the 
radiative force due to gravity darkening is calculated and 
in §4 the drift velocity of metals around an equipotential 
is derived. Numerical results of the diffusion velocity in 
stellar envelopes are presented in §5, and by considering 
the relevant timescales, §6 identifies the important regions 
within a star for diffusion to occur. A discussion is pre- 
sented in §7, and conclusions given in §8. 
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2. Transport of chemicals 

The conservation equation for a trace element of concen- 
tration c is 



p-^ + V. {-pD. Vc + puc + pvc} 



(1) 



where D is the diffusivity tensor, u is the global motion 
of the fluid (in this case the meridional motions gener- 
ated via the rotation of the star), and v describes any 
motion affecting only the contaminant. The concentration 
c is defined as the number density of ions compared to the 
number density of protons. Hereafter, the star is assumed 
to be axisymmetric, so that using a spherical coordinate 
system, for any X, dX/d(f> = 0. Also, for rotating stars, the 
cquipotcntial is not defined by r —constant. Therefore the 
"horizontal" direction is taken to be perpendicular to the 
equipotential, and the "vertical" direction to be parallel 
to equipotential. 

Following the work of Chaboyer & Zahn (1992) and 
Zahn (1992), Charbonncau (1992) examined the effect 
of anisotropic turbulent diffusion on the distribution of 
c(r,0) (Zahn 1975, Tassoul & Tassoul 1983, Zahn 1987), 
and found that if the horizontal Reynold number is Rr < 
1 then the equation above loses its bi-dimcnsional be- 
haviour, and becomes one dimensional ie. c(r,6) — ► c(r). 
Essentially the horizontal turbulence becomes vigorous 
enough that it may redistribute the concentration of the 
trace element faster than meridional motions can generate 
an anisotropy. This also makes the angular momentum in- 
dependent of 9 and gives rise to "shellular" rotation of the 
star (Zahn 1992). 

In the current study, however, the velocity field v may 
be non-zero due to the meridional forces on trace elements 
arising from the anisotropy of the radiation flux (see §2). 
To assess its likely effect in the presence of bulk meridional 
motions and diffusion, the ratio of the (turbulent) diffu- 
sive term to the radiation driven advective term across an 
equipotential s in cq.l is examined: 
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where Dh is the horizontal diffusion coefficient, the extra 
velocity component is v = (0,i> s ,0) and h(c) is a func- 
tion containing logarithmic derivatives of the concentra- 
tion with respect to 9. This is only an approximation as 
the unit vectors along an equipotential s and 9 are not 
parallel. However, this expression will suffice for the cur- 
rent study. Also it should be noted that there is no term 
in the angular momentum evolution equation equivalent 
to the radiation term V(pvc) in eq.l. Therefore, shellu- 
lar rotation will still be achieved in the star, although the 
distribution of the trace element may be anisotropic. 



The discussion of turbulent diffusion coefficients Dh 
has a long history, and the simple approximation due to 
Chaboyer & Zahn (1992) is adopted here: 



D H = 



\ru r \ 
C H 



(3) 



where Ch is a number of order unity and u r is the radial 
component of the velocity of the bulk meridional circula- 
tion. With this replacement eq.2 becomes 
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v s C H 



h(c) 



(4) 



The main controlling parameter here is the ratio of the 
radial component u r of the bulk meridional flow and the 
velocity v s of the trace element due to the radiation as it 
slips through the plasma. 

If 1Z > 1 then the turbulent diffusion term will domi- 
nate the movement of trace elements, and so the star will 
become chemically homogeneous on equipotentials. How- 
ever, if 72. < 1 then the distribution of the elements is 
determined by the extra drift velocity vg. The criterion 
adopted in this paper for chemically inhomogeneous stars 
is that the extra drift velocity v s is larger than the typical 
velocity of meridional motions u r . In a recent numerical 
study of the a 2OM star, Urpin et al. (1996) found that 
for high rotational velocities, the circulation velocity is less 
than 3 x 10 _5 cm s^ 1 . This value is significantly lower than 
that from Eddington-Sweet theory, which fails close to the 
surface of the star. If velocities of v s > 10" 5 cm s _1 are 
produced due to the effects of rotation, then the star may 
develop metal-rich and metal-poor regions on the same 
equipotential. 

3. Radiation forces 

3.1. Flux around the star 

Assuming that the potential of the outer parts of a rotat- 
ing star may be represented using the Roche approxima- 
tion, the radius of equipotential surfaces r is dependent 
on the angle to the pole 9, and is found by the solution of 
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where ui is the angular velocity, r po i is the radius in the 
polar direction, and M is the mass contained within the 
equipotential surface. This has been solved for r(9) geo- 
metrically by Collins & Harrington (1966). Deep within 
the star the Roche approximation will be incorrect, al- 
though for the present discussion, it is adequate. Indeed 
the region in which it is possible to generate latitudinal 
metallic abundance anisotropics is found, a posteriori, to 
be the outer parts of the envelope and soM« M* , and 
eq.5 will be a good representation. 



John M. Porter: On the possibility that rotation causes latitudinal abundance variations in stars 



3 



The flux of the star / is proportional to the local grav- 
ity (von Zeipel 1924): 



/ocV 



GM _ cu 2 r(0) 2 sm 2 6 
H0) 2 



(6) 



Along the surface of any equipotential the local flux / is 
higher at the pole than the equator (gravity darkening). 
This difference between pole and equator is the root cause 
of the latitudinal variation of the metallic abundance. It 
is noted that this gravity darkening actually causes the 
meridional circulation, else energy would not be conserved 
in the star. However, the change in flux is not destroyed 
by the motions. 

3. 2. Radiation force on a trace element 

The force on ions due to radiation g\^ A is dependant on 
the gradient of the flux passing through that point, i.e. 
g* ad oc V/. To facilitate the calculation of the force along 
an equipotential the ratio of the forces in the radial and 
latitudinal directions is now calculated. In the outer parts 
of the star, the luminosity is a constant, and hence the 
flux varies with distance r from the centre of the star as 
/ oc r~ 2 . As the radiative force is proportional to V/, 
then the ratio of the latitudinal to radial components is 



9lf 



r din/ 

2 ds 



(7) 



where s is the distance along an equipotential. Again, note 
this is only an approximation as the unit vectors s and f 
are not orthogonal. This ratio is shown in fig.l for different 
rotation rates. As can be seen it is typically 0.1 — ► 1, and 
therefore, perhaps surprisingly, the radiative force around 
an equipotential is typically similar to the radial radiative 
acceleration. 

Michaud et al. (1976) see also Alecian & Artru (1990) 
have discussed the radiative force due to lines on heavy 
elements. Their expressions, along with eq.7, lead to the 
latitudinal force on on a trace clement (with a negligible 
gradient in the concentration of the ions) of 



-5.574x 10 



7 1 e 4 

A,T 4 



R„ 



r din/ 

2 ds 



(8) 



where T4 and T e 4 are the temperature and effective tem- 
perature in units of 10 K respectively, r is the radius and 
i?* is the stellar radius. The atomic parameters are en- 
veloped in <&i, and Ai is the atomic mass of the element 
i. 

With the inclusion of a non-zero concentration of ions 
the driving lines become saturated, and the line force is 
reduced from eq.8: 
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Fig. 1. Ratio of latitudinal to radial radiative force defined 
from eq.4. The lines correspond to a rotation of 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8 and 0.9 of break-up velocity, where the uppermost 
line corresponds to the highest rotation, decreasing downwards. 



(Alecian 1985). The reference concentration cq is 

c = 9.83 x 10~ 23 fin e T~ 1/2 y 2 (10) 

where R is the Rosseland mean opacity, n e is the num- 
ber density of electrons and ^ is an atomic parameter 
containing line broadening effects. 

4. Drift velocity 

To calculate a typical order of magnitude of the drift ve- 
locity v between the ions and rest of the plasma, the radia- 
tion force will be approximately balanced by the retarding 
frictional force. The diffusion velocity is 



D12 -Vine + A., 



kT 



rad 



(l-2A t + Z t )^-a T V\nT) 



(11) 



Here, m p is the mass of a proton, k is Boltzmann's con- 
stant, Zi is the charge of the element i in units of the fun- 
damental charge e. The thermal diffusion coefficient is 
and the microscopic diffusion coefficient is Z? 12 . The effec- 
tive gravity is zero in the s direction (from the definition 
of the equipotential), and so the third term in parentheses 
above is zero. Also, the thermal diffusion term is typically 
small compared to the radiative term and so the diffusion 
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velocity is determined only by the radiative forces, and 
the concentration gradient. 

If the mass of species 1 (the stellar plasma - mostly 
hydrogen and helium) is neglected compared to the mass 
of species 2 (the diffusing ions), and ion shielding is ne- 
glected in the computation of the Debye length, eq.40 of 
Paquette et al. (1986) becomes 
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where p is the density. To obtain a simple expression for 
the diffusion velocity, an initially chemically homogeneous 
star is assumed, so that Vc = and eq.9 are good repre- 
sentations. The diffusion velocity is then 



-2.2 x 10" 
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(14) 



This expression will overestimate the drift velocity as soon 
as a concentration gradient is established, as then Vc ^ 
in eq.ll. However, with these points in mind, eq.14 is used 
as an estimate for the velocity in eq.4. 

5. Numerical estimates of the diffusion velocity 

Here the typical drift velocity is calculated for two model 
stars. The stellar models have been provided by M. Salaris 
(private communication) and are coomputed according to 
the evolution code described in Salaris et al. (1997) and 
references therein. The two models correspond to (2.0M Q , 
l.8R Q T cS = 8700K) and (5.OM , 3.Oi? T cff = 16500K) 
main sequence stars. These models are computed with no 
rotation and solar composition. Although below the mod- 
els are used to represent non-solar composition rotating 
stars, it is unlikely that this slight inconsistency will in- 
troduce a large error. 

The calculation of the atomic parameters and ^ 
are lengthy and cumbersome, and so the parameters given 
by Alecian et al. (1993) for iron are assumed (their ta- 
ble 4), derived using Opacity Project data. These are 
for the states FelX- FeXVII. For a given temperature in 
the model stars, the atomic parameters used are inter- 
polated between those for the dominant, and next domi- 
nant ions (the temperatures at which ions are dominant is 
given in eq.3 of Alecian et al. 1993). These temperatures 
also give the applicable range used in the calculation : 
1.95 x 10 4 -ftT < T < 6.0 x 10 5 K to cover all the ions used. 

Although the models are calculated with solar metal- 
licity, the drift velocity v s is calculated for two values of 
iron abundance: solar abundance c = 3 x 10~ 5 and 10 -2 



solar. Although strictly this is inconsistent, the error as- 
sociated in doing so will be very small. 

In fig. 2, the radial radiative acceleration is shown for 
the two stars with solar iron abundance, and should be 
compared with fig. 4 of Alecian et al. (1993). It can be seen 
that the interpolation used here to calculate the atomic 
parameters is not introducing any large errors in the cal- 
culation. 
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Fig. 2. Numerical estimates of the radial radiative acceleration 
g[ a r d for the two model stars. 



Fig. 3 shows the value of v s as a function of the frac- 
tional mass exterior to that point AM/M for the two mod- 
els. The stars is assumed to rotate at 0.7 of their break up 
velocities (0.89 of its break-up angular velocities), equal 
to v Iot = 260 km s _1 and 320 km s _1 respectively for the 
2M Q and 5M Q stars. The top panel of fig. 3 corresponds 
to the 5M@ star, and the lower to the 2M star. In each 
panel, the lines refer to 9 = 10°, 30°, 50° and 70° with 
the uppermost line corresponding to 70° decreasing down- 
wards. The solid and dotted lines correspond to iron abun- 
dance of solar and 10 _2 solar respectively. 

In three of the four models, the drift velocity exceeds 
10~ 5 cm s" 1 in the outer layers at some latitude. This in- 
dicates that there may be regions in the envelope where 
iron may start to diffuse round the star. For the 5M star 
the outer 10" 5 ' 5 (10~ 4 ' 7 ) of the solar (10~ 2 solar) com- 
position model may suffer latitudinal drift of iron. For 
the 2M star, at 10~ 2 solar composition the situation is 
similar - w s > 10 _5 cm s _1 in the outer 10~ 5 of the star. 
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Fig. 3. Numerical estimates of v 3 from eq.14. The four lines 
in each panel refer to velocities at 6 — 10°, 30°, 50° and 70° 
with the uppermost line corresponding to 70° decreasing down- 
wards. The solid (dashed) lines correspond to solar (1CF 2 solar) 
abundance of iron. 



However for solar composition, there may be no effect as 
the radiative acceleration on iron is smaller. In this case 
the turbulence may homogenise the star faster than the 
metals may drift due to the radiation. Taken at face value, 
fig. 3 indicates that stars rotating at significant fractions of 
their break-up velocity may be prone to latitudinal metal- 
lic drift in their outer envelopes. 



6. Timescales 

Is is likely that a star may come into equilibrium? This 
clearly depends on the timescales over which the relevant 
processes act. Here the important timescales are (i) the 
diffusion time (time taken for the ions to diffuse round 
that star), (ii) the main-sequence lifetime of the star, (hi) 
the lifetime of each layer (if mass-loss is present). 

6.1. Diffusion, stellar and mass-loss timescales 

For an asymmetry to occur in the elemental distribution 
along an equipotential, the metals must diffuse a signif- 
icant distance from their initial latitude to the equator. 
This therefore defines a drift timescale td, dependent on 
latitude it/ 2 — 9 , which is that over which a latitudinal 
asymmetry in the metals will be generated: 



td{6) 



R« (tt/2 - 



(15) 



The star only has a finite time in which to allow met- 
als diffuse around its equipotentials. The main sequence 
lifetimes for these models are tms = 9.5 x 10 8 yr and 
7.9 x 10 7 yr for the 2M Q and 5M Q star respectively. 

The mass-loss timescale within the envelope depends 
on the mass exterior to that point (AM) and the mass 
loss rate M. If the mass-loss rate M is a function of 
azimuthal position 9 for rotating stars (for line-driven 
winds see Friend & Abbott 1982 and for radiation-driven 
dusty winds see Dorfi & Hofner 1996), then the mass-loss 
timescale also varies with 9. However, as a first approach 
tm{9) is taken to be independent of 9: 



tm 



AM 
M ' 



(16) 



6.2. Where can the asymmetry exist? 



In order for any significant metallic asymmetry to be gen- 
erated, the diffusion timescale must be less than the mass- 
loss timescale. This ensures that the diffusion process has 
enough time to act before that layer of the star is lost in 
an outflow. The diffusion timescale must also be shorter 
than the typical main-sequence timescale else no signifi- 
cant asymmetry will build up over the lifetime of the star. 
Therefore, chemical abundance variations may be gener- 
ated in regions for which the two inequalities 



TD^T M , T D (6) £TMS 



(17) 



are satisfied. 

The timescales are displayed in fig. 4. The lines are 
the diffusion timescale around the star for 9 = 10°, 30°, 
50° and 70° with the lower line corresponding to 70° de- 
creasing upwards. The solid lines are for solar iron abun- 
dance, and the dotted lines are for 1CP 2 solar abundance. 
The oblique dashed lines are the mass loss timescale, cor- 
responding to mass-loss rates of 10 _12 M Q yr _1 (lower), 



6 



John M. Porter: On the possibility that rotation causes latitudinal abundance variations in stars 




-6 -5.5 
log 10 ( AM / M ) 




4.5 



( AM / M ) 



Fig. 4. The timescales in the star: the diffusion timescale 
around the star for 6 = 10°, 30°, 50° and 70° with the 
lower line corresponding to 70° decreasing upwards. The 
solid (dashed) lines correspond to solar (1CP 2 solar) abun- 
dance of iron. The dotted lines are the mass-loss timescales 
corresponding to mass-loss rates of 1CF 12 Mq yr _1 (lower), 
1O- 13 M yr" 1 , KT^Moyr- 1 (upper). 



l(r 13 M© yr" 1 , l(r 14 M yr" 1 (upper). Fig.4 should be 
interpreted as follows: the region to the right and lower 
than the mass-loss timescale lines may, if the diffusion 
velocity is larger than the diffusive velocity (~ the merid- 
ional velocity u r see §2) develop metallic abundance gradi- 
ents around an equipotcntial. Therefore, the 5M model 
may have enough time to develop a significant abundance 
gradient across most of the meridian in the outer regions 
is the mass- loss rate is ~ fCT 13 M Q yr -1 . 

Fig.4 suggests that during the main sequence lives of 
the 10 _2 solar abundance 2Mq and both 5M© stars, the 
drift velocity due to the anisotropic flux in the star v s is 
large enough that a latitudinal abundance gradient may 
be generated. However, this statement must be confined to 
stars which are rotating at significant fractions of break- 
up, and only applies to bands near to the equator. 

Although this effect may produce metallic drift, it 
seems from fig.4 unlikely that the whole of an equipoten- 
tial s will be involved. It appears that only regions close to 
the equator will be effected and consequently the metal- 
lic distribution will not come into equilibrium during the 
star's main-sequence lifetime. The distribution of metals 
is therefore difficult to calculate - indeed the dynamical 
modelling of the distribution is out of the scope of this 
paper and is flagged for further study. 

7. Discussion 

The results of the previous sections arc startling. It has 
been found that in bands near the equator metals may 
drift toward the equatorial plane due to the anisotropic 
radiation field created when the star rotates. However, it is 
noted that if the mass-loss rate in the equatorial regions is 
strongly dependent on the metallicity, then the enhanced 
mass-loss rate will reduce the time that a given layer is 
bound to the star. Therefore the metallicity enhancement 
will not be able to build up to such a large level (as esti- 
mated below). In this case then a new equilibrium will be 
approached in which the magnitude of the metallicity en- 
hancement is controlled by the enhanced mass-loss. This 
aspect of the problem is currently being investigated - the 
homogenising effect of the outer convection zone in the 
models has, so far, been neglected. This may homogenise 
the metallicity fast enough that no effect on the mass-loss 
rate may be observed. 

The drift will slow and eventually stop when the loga- 
rithmic concentration gradient becomes comparable with 
the radiative force (see eq.ll): 



dine 



1 kT 



(18) 



(again note this is an approximation as s and 6 are not par- 
allel). As an example, the temperature is set to 2.5 x 1 5 K 
corresponding to AM/M ~ f CP 6 for both stars. The radii 
of the stars of r = 3.C1R© (5M ) and 1.8-R© (2M Q ). 
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From fig. 2 the radial radiative force is w 10 4,7 (10 4 ' 1 ) for 
the 5M© (2Mq ) star. If it is assumed that the radiative 
force around the star <?™ d is a tenth of this (see fig.l) 
then the right-hand side of eq.18 is 1360 and 205 for the 
5M Q and 2M© stars respectively. These are clearly large 
numbers - indicative of a very large abundance build-up in 
the equatorial plane. However as shown in §6 it is unlikely 
that the whole of the meridian will come into equilibrium. 
It may be surmised, though, that significant abundance in- 
homogeneities can be generated (without being thwarted 
by the concentration gradient) through rotation. 

In the course of this paper, radial motions due to ra- 
diative acceleration or gravitational settling have been ne- 
glected. It is likely that a radial component of drift velocity 
is present as well as the meridional component focussed 
on here. Therefore, the motion of heavy element ions will 
not be solely be meridional - if the radial radiative accel- 
eration exceedes gravity, then the drift will be towrd the 
equator and toward the surface. 

It is difficult to accurately assess the effects of large 
latitudinal abundance gradients on the structure of a star. 
As soon as the abundances change at a given point, then 
the flux distribution will also change along with the local 
convective stability criterion. Convection smooths out the 
abundance overdensity, and provides some feedback to the 
ionic build up. It is posssiblc that this feedback will regu- 
late the ionic drift, modifying the structure of parts of the 
star. 

Let us now consider the evolution of the outer enve- 
lope of the star during the main sequence. As the very 
outer layers are lost in a wind, then a given layer be- 
comes slightly more diffuse and the layer moves to lower 
AM/M (which will promote latitudinal ion drift). How- 
ever, as the layer becomes convectively unstable then 
it is quickly made homogeneous. The presence of extra 
radiation-blocking ions in fact make the layer more un- 
stable to convection, and so the equatorial regions of ro- 
tating stars can be expected to have a more extensive 
convective region. Due to this rapid convective homogeni- 
sation, it is very unlikely then that this latitudinal drift 
will be observed in the photosphere of the star during its 
main sequence lifetime (the possibility of a star having 
different chemical compositions on the same equipoten- 
tial has already been considered using radial diffusion and 
magnetic fields by Michaud, Megessier & Charland 1981). 
The possibility of observing a direct manifestation of lati- 
tidinal abundance variation in post main-sequence phases 
of evolution is currently under study. Here the mass-loss 
rate may be large enough to prevent complete convective 
chemical homogenisation of regions which have generated 
abundance gradients during the main sequence. 

It is clear that as the timescales over which a significant 
abundance gradient may be generated in parts of the star 
are similar to cither mass-loss timescales or main-sequence 
lifetimes. Therefore the evolution of the star needs to be 
taken into account whilst considering the abundance time 
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dependence - something which is not attempted here. Al- 
though difficult due to necessarily 2D nature of the calcu- 
lation, it appears that time dependence of the abundance 
distribution this needs to be included in models for rotat- 
ing stars. 

8. Conclusion 

It has been shown that rotation may indirectly have a 
large effect on the latitudinal abundance distribution of 
metals via gravity darkening. Indeed, the effect is possibly 
active in all the model stars considered here. Stars of low 
metallicity are particularly prone to this effect. If the mass 
loss of the star is strongly metallicity dependent, then the 
effect in the outer layers may be curtailed somewhat. This 
may lead to an extra mass-loss rate enhancement in the 
equatorial plane. 

Although several major points regarding the effects of 
meridional diffusion of ions toward the equator have sim- 
ply been touched on, the principle of metallic latitudinal 
distribution asymmetries has at least been put on quanti- 
tative basis. 
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